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Abstract: The paper describes the design and 
experimental characterisation of a small 
reciprocating linear permanent-magnet generator 
which has been developed for on-board 
generation of electrical power for telemetry 
vibration monitoring systems. A comprehensive 
design methodology, which embraces both the 
electromagnetic and mechanical design, is 
described. The electromagnetic design synthesis is 
based around analytical field solutions which 
allow the optimisation of the generator 
parameters, and is further verified by nonlinear 
electromagnetic finite-element analysis. Nonlinear 
structural finite-clement analysis is employed to 
aid the design of flexural disc springs with the 
desired stiffness and stroke. The utility of the 
design methodology is illustrated by a design 
study and extensive experimental charact- 
erisation of a prototype device. 



List of principal symbols 

A = magnetic vector potential, Wb/m 

B = magnetic flux density, T 

Brem = remanencc, T 

/ = frequency, Hz 

H = magnetic field strength, A/m 

/ = output current, A 

J = winding current density, A/m 2 

K E = EMF constant, V-s/m 

K T = force constant, N/A 

K, = spring constant, N/m 

L - winding inductance, H 

M = remanent magnetisation, AM 

m - mass of plunger, kg 

P„ = output power, W 

R L = load resistance, Q 
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plunger radius, m 
stator radius, m 
axial displacement, m 
permeability of free space, H/m 
relative recoil permeability 
pole pitch length, m 
angular frequency, rad/s 
winding flux linkage, Wb 



1 Introduction 

There are an increasing number of applications which 
demand the remote provision of low levels of electrical 
power and for which on-board power generation is 
preferable to the use of batteries, which have a limited 
lifetime and contain toxic materials [1]. For example, in 
telemetry-based condition vibration monitoring sys- 
tems, where there is potential to utilise the inherent 
vibrations to generate the required electrical power via 
a reciprocating linear generator. In principle, any of the 
linear reciprocating actuator topologies that have been 
developed for use in applications ranging from auto- 
motive [2] to healthcare [3] may also be used for power 
generation. However, tubular machines with perma- 
nent-magnet excitation have a number of distinctive 
features [4], such as a high power density and a high 
efficiency, which make them particularly attractive for 
reciprocating linear generators. 

This paper describes the design and experimental 
characterisation of such a small reciprocating linear 
permanent-magnet generator, shown in Fig. {. It is a 
single-phase tubular device, which is capable of pro- 
ducing 20mW at 5V, with a nominal resonant fre- 
quency of 50H7. and a stroke of ±0.8mm. The two 
axially-magnetised sintered NdFeB magnets and the 
mild-steel pole pieces give rise to an essentially radial 
magnetic field in the region occupied by the generator 
winding. The moving permanent-magnet plunger is 
supported at each end by beryllium copper disc springs, 
which provide a high degree of radial stiffness, and an 
axial stiffness that can be accurately controlled by 
appropriate design of the spiral grooves. This spring 
plunger assembly mechanism also eliminates mechani- 

-xai- wc ar, a nd there fbxe-xesult.s i n i mp xo^d-jreliahility 

[5]. When the generator is vibrated, as shown schemati- 
cally in Fig. 2, electrical power is generated as the 
result of the relative movement between the plunger 
and the winding. For a given amplitude of vibration, 
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the system output power reaches a maximum when the 
frequency of the vibration matches the mechanical res- 
onant frequency of the system. 




o 

Fig. 1 Reciprocating tulmlar linear permanent-magnei generator 
a Linear generator " 
b Fkxural spring 



windings 



permanent 
magnat-plunger 




Fig. 2 Schematic of reciprocating linear generator system 

To facilitate design optimisation and accurate 
dynamic modelling, a variety of techniques have been 
employed to predict the magnetic field distribution in 
tubular linear permanent-magnet machines [6], the 
most common approach being to employ a lumped 
equivalent circuit [7, 8J. While this allows the relation- 
ship between critical design parameters and machine 
performance to be established analytically, it suffers 
from problems associated with model inaccuracy, par- 
ticularly when flux leakage is significant and the flux 
paths are complex. Therefore, numerical analysis of the 
field distribution and evaluation of performance [9, 10] 
is also employed. However, while techniques such as 
finite-clement analysis provide an accurate means of 
determining the field distribution, with due account of 
saturation etc., they remain time consuming and do not 
provide as much insight, as analytical solutions into the 
influence of the design parameters on the machine 
behaviour. 

To overcome these problems, an analytical solution 
for the magnetic field distribution has been established 
in the cylindrical co-ordinate system [1 1]. It is an exten- 
sion to the analysis which has been developed for tubu- 
lar linear permanent magnet machines to account 
specifically for the relatively short stroke of reciprocat- 
ing linear machines. The analytical solution allows th e 



and system dynamic modelling. The predicted results 
are verified by finite-element analyses and measure- 
ments on a prototype generator. 

2 Field distribution due to permanent-magnet 
source 

2. 7 Magnetic field distribution 
Fig. 3a shows the longitudinal section of the linear 
machine in which the axial length of the iron sleeve is 
usually longer than that of the plunger. To establish an 
analytical solution for the magnetic field distribution 
the following assumptions are made: 

(i) The axial length of the machine is infinite and com- 
prises an infinitely long iron sleeve and a series of seg- 
mented plungers extending to infinity along the 
machine axis, as shown in Fig. 3b. 

(ii) The armature is slotless, and the permeability of the 
iron sleeve and pole pieces is infinite. 

("0 The relative recoil permeability of the magnets u,. is 



6f* 




J 3EDHT 

region til 



u 

Fig. 3 Analytical field models 

a Field model 

* Extended field model 

The magnetic field model of the actual linear genera- 
tor may be represented by one of the repetitive ele- 
ments in Fig. 3b, provided that the separation distance 
t,. between two adjacent plungers is much greater than 
the pole pitch x p , i.e. r r » r p . Consequently, the mag- 
netic field analysis is confined to three regions, viz. the 
airspace and winding regions J and III in which the 
permeability is ^, and the permanent-magnet region // 
in which the permeability is HqUt- To simplify the anal- 
ysis, the regions II and III are combined by assuming 
that fx u = n m = lifj, i.e. the airspaces between the 
plungers are filled by permanent-magnet material, but 
the magnets are only fully magnetised in region II, 
being unmagnclised in region III. Therefore 



B = 



in region / 



1 Mo/J-H + MoM in regions 77 and III 



prediction of the thrust force and back-EMF in closed 
forms, which facilitates the characterisation of 
machines and provides a basis for design optimisation 



(1) 

where M is the remanent magnetisation. For a perma- 
nent magnet having a linear demagnetisation character- 
istic, n r is constant and the magnetisation Ml is related 
to the remanence, B rm by M = IW/4>- Further, for a 
rare-earth magnet p r s 1.10. Hence, the assumption 
lhaL-M,. = 1 0 is justified . It is convenient to fo rmulate 
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the field distribution in terms of a magnetic vector 
potential A defined as B = V x A, and the cylindrical 
co-ordinate system shown in Fig. 4. The governing field 
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equations, in terms of the Coulomb gauge, V A = 0, 
are 

f V 2 A/ =0 in region 7 

I V 2 A//jn = -MoVxM in regions 77 and 777 

(2) 




Pig . 4 Cylindrical co-ordinate system 



Since the field is axially-symmetric, A only has the 
component Ag, which is independent of 0. This leads to 

in region 7 

= ~^oV x M 
in regions 77 and 777 
(3) 

In "the cylindrical co-ordinate system, the flux density 
components arc deduced from A e by 



£ (\ 8 
dz 



The magnetisation M of the linear machine is 
M = Af,e. 



Wmo 



(4) 
(5) 



Fig. 5 Magnetisation distribution 



(t p -i m )/2 (tp+T m V2 



where M z as shown in Fig. 5 denotes the component of 
M in the z-direction, and may be expressed as the Fou- 
rier series 



Mz = V !^ S i n taA 



?17TT„ 

sm I — 1 sin m n z 



(6) 

where x m is the axial length of the magnet, x lp = x p + xj 
2, and m„ = wJt/r^. The boundary conditions to be sat- 
isfied by the solution to eqn. 3 are 

. Bj.C^ =0; B rT . ni r L_ n = 0; 

S//r| N=(r)>±Tm)/2 = 0; J?//;.| w=(rii) _ T ,. /2) =0 

Brzl^Rr,. = 0; 5 ; Jr=fl m = o 



B j r \ r—Rw 

' '(T„-T m )/2SI*l<(T, + T m }/2 



77/Jr=« m 

'(t p -t ^)/2<| 2 |<(T,H-T, n )/2 



('■p-r„,)/!<kl<{r,+T„,)/2 



= 7/ /; Jr=fi„, 

1 (-,.-r m )/2<|^|<(r p+Tm )/2 



'»><M<T,p '^<|S|<T| F 

f 

Jo 



H !z \r=n,„ =Huj z \r=n.„ 



=(r p -T m )/2- 

(T>-T m )/2 



Jo 



(7) 

Solving eqn. 3 subject to the boundary conditions of 
equ. 7 yields 

B lT (r,z) 

CO 

= ~ [a/n5Ji(m n r) +6 J „S7i' 1 (m n r)]cos(m„z) 

n=l,2,... 

DO 

= l a ^ BJ o( Jn nr) - bj n BK Q (m n r)] sm(m n z) 

n=l,2,... 

(8a) 

Bi Ir (r,z) 

CO 

i=l,2,... 



(86) 



= l a "'- 57 o(?//» r )] cosq;n(z - m) + Bo 

>=»1,2,... 

Birir(r,z) 

CO 

= X) l a /W 5J 'i(9///i I ')]sm«'7//i(2 - t p ) 
j=r,2,... 

Bnn(r, z) 

CO 

(8c) 

where B1q('), BI\{-) are modified Bessel functions of the 
first kind; BK 0 ('), S7T|(') are modified Bessel functions 
of the second kind, of order 0 and 1, respectively; q IH , 
q mj , and x u are given by 

qui = m/r m ; m = (r p - r m )/2 

<•/;;/., = (2j - l)it/T T 

and a,„, £>/,„ <7 /Wy -, and B 0 are defined in the Appen- 
dix 

~22 ~ Comparison with finite-element" " 
calculations 

The main design parameters of the reciprocating linear 
machine for which analytical field solutions have been 
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obtained are given in Table 1. The magnets are sintered 
NdFeB with B r ,,„ = 1.2T. The analytical field distribu- 
tions have been validated by finite-element calculations 
of the radial and axial variations of the radial and axial 
flux density components in the airgap/winding 
regions. " 

Table 1: Parameters of reciprocating linear permanent- 
magnet generator 



R, (mm) 


8.5 


fii(mm) 


7.0 


R„ (mm) 


6.0 


t„, (mm) 


11.0 


t p (mm) 


11.0 


*W (mm) 


7.5 


r m (mm) 


7.5 


r,(mm| 


44.0 



The finite-element solutions were obtained by apply- 
ing a periodic boundary condition at the axial bounda- 
ries z = ±r, p and imposing the natural Neumann 
boundary condition at the surfaces of the stator sleeve 
and plunger iron pole pieces. Fig. 6 compares flux den- 
sity components as functions of axial position z al a 
constant radius r = 7.0mm. It will be seen that the ana- 
lytical solutions agree extremely well with the finite-ele- 
ment results. 




0.01 0.02 0,03 
axial position z,m 
a 




J.-O.0S \ - 



0.01 0.02 0.03 
axial posMon z,m 
b 

0 007m W " X dav ^ y com P° ne "' 1 as functions of axial position at r = 
a Radial component 
b Axial component 
— — analytical 
— — FEM 

3 EMF and force prediction 

3. 7 Flux-linkage and EMF 
Assuming that each coil of the armature winding com- 
prises a number of circular turns, and occupies an area 
bounded by ;• = R h r = R„ z = z - and z 2 = z + 
rJ2, as shown in Fig. 7a, where t„ is the coil axial 
width, the coil flux-linkage may be obtained from the 
following integration: 

ft, 



= J J 2irrA ro (r,z)drdz 



J-T,„/2 it, 



/, *n sin m n z 

o=l,2,... 



(10) 



where 

$„=^sin(m n ^) 
ft.. 

x J r [ajnBIj (m»r) + bi n BK, (m n r)] dr 

(ID 

The total flux-linkage of the distributed armature wind- 
ing, comprising coils 1, 2a and 2b, as shown in Fig. lb, 
is therefore given by: 

CO 

i>wp — tpwi - (V w2a + Tpw2b) = *n P sin m n z 

n=l,2,... 

(12) 

























1 

•' ft 










1"' 
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Fig. 7 Distribution of armature winding 
a Cotl distribution 
b Winding distribution 

where 

4> = ^ n K ipn N wv 

ft. 

x J r [a ln Bh (m n r) + b In BKi (m n r)) dr 
fti 

(13) 

K dp „, defined as the winding factor of the nth harmonic, 
is given by 



K, 



sin(m n r Ml /2)-2sin(m n r v2 /2)coB[m, l (r wl +T v , 2 )/2} 

(14) 

and N wp is the number of series turns of the winding. 
Hence, the winding induced EMF is obtained as 



„ _ r H'w P 



~ ^2 KiSn. cos m n z ) ~ 

n=l,2,... ' 



(15) 

where K En is analogously defined as the back-EMF 
constant of the nth harmonic, and is 
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R, 

x I r {ai n DI x (m n r) + b ln BK x (m n r)} dr 
k, 

(16) 

3.2 Thrust force 

With reference to Fig. 7a, the thrust force exerted on 
the coil when it carries a current density J may be 
obtained from the following integration: 

i + T w /2 R, 

?n, = - J I ' 2nrJ£, r (r,z)drdz 
which may be written as 



where F„ is 



F ™ = Yl FnCQsm n z (17) 

n=i,2,... 



m n V 2 J 
R, 



* J r [a ln Bh Kr) + (m n r)] rfr 

(18) 

The total force exerted on the distributed winding a 
carrying current I is therefore 

Fwp - ^ A> n cos m n ^i (19) 

where ^ is analogously defined as the force constant 
of the rcth harmonic and is given by 

(R. - Ri) 

R, 

x J r{a In BI l {m u r)+bi n BK l (m. n r))dr 

Ri 

(20) 

which is identical to the EMF constant of the nth har- 
monic K E „. 



0.06 




axial plunger position z, mm 

F'S- 8 an f^ t f c '/f' ! * ast ' ""/""Clion of tixial plunger displacement 

— O — measured 
— *— FEM 
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3.3 Validation 

Fig. 8 compares the measured winding flux-linkage of 
the single-phase reciprocating linear generator shown in 
Fig. 1«, whose design parameters were given earlier in 
Table 1, with the analytical and finite-element calcula- 
tions. Both saturation and end-effects associated with 
the finite length of the stator sleeve are now accounted 
for in the finite-element calculation. However, as will 
be seen, the analytical prediction agrees well with both 
the finite-element and measured results, thus validating 
the developed analytical model for the linear machines, 
for which the computational times are significantly 
shorter than those required for corresponding finite-ele- 
ment analyses. 

4 Design optimisation 

4. 7 Frequency response 

As seen in Fig. 8, if the deflection z is small such that 
sinKz) = m„z, (n = 1, 2, AO, where N is the 
number of harmonic terms used for the calculation of 
the flux-linkage, y wp is linear with respect to z. Thus, 
cqns. 15 and 19 may be simplified as 

dz 

c-w 7 , = —K E -r-\ 



'dt' 



F wp — Kti* 



(21) 



(23) 



where K E and K T are defined as the EMF and force 
constants of the machine, respectively, and arc 

N 

K B = K T = I<Bn (22) 

n=:l,2,... 

Fig. 2 shows a simplified mechanical model of the lin- 
ear power generation system. Assuming that the mass 
of the vibrating member, to which the generator bous- 
ing is attached, is much greater than the mass of the 
plunger, so that the source of vibration is not affected 
by the plunger movement, the governing equations of 
the electromechanical system may be written as 

Lf t +(R + R L )i = -K B f t 

where m is the plunger mass, K, is the stiffness of the 
suspension spring, K D is the damping factor, L and R 
are the inductance and resistance of the generator, 
respectively, and R, is the load resistance. For a sinu- 
soidal excitation vibration x = V2X sin wr the steady- 
state plunger deflection Z and output current / as func- 
tions of the vibration angular frequency to are given, 
respectively, by 

zh = ujWm 

where R T = R + R L . Fig. 9 shows Z, / and the output 
power P„ = R L f as functions of the vibration fre- 
quency/, assuming X - 0.13 mm, m = 0.020kg K F - 
66.7 V-s/m, L = 0.15H, R = I.5kG, R L = 1.5kQ, K D = 
0, K ? = 2000N/m. The plunger deflection Z reaches its 
maximum amplitude approximately at the resonant fre- 
quency co c = V(KJm), which corresponds to the maxi- 
^jjmjqiUpuLrjica vr i r for a giv e n vibration a mplitud e, 

4.2 Winding design 

Since the generator employs a slotless armature the 
winding has a very low LIR ratio. Therefore eqn. 25 
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may be simplified as 



R + Ri 



Z 



(26) 



< 

c" 




frequency, Hz 

Fig. 9 Plunger deflection and output as functions of vibration frequency 
a Plunger deflection 
b Output current 
c Output power 

From eqns. 16 and 22, the EMF constant K E is propor- 
tional to the number of winding turns and may be 
expressed as 



where 

N 



Kb = K W N WV 
2irK d 



(27) 



Mp 71. 



(R, - Ri) 

J r [ajnBh (m n r) + b ]n BKy (m n r)} dr 



n=l,2,... 
ft. 



(28) 

The winding resistance R is also related to N wp as fol- 
lows: 



where 



A' r = 



R=K T h% 
w(R, + Ri)p c 



Pf(R s -R l ){r wl +2r w2 ) 



(29) 
(30) 



p cu denotes the resistivity of copper, and P f the winding 
.jacking factor. Thus, th e output power to a resistive 
load R L may be expressed "as " " 



P„~I 2 Rl 



Ky,N wp Zu 



K r Nl p + R L \ 



Ri (31) 



It can be shown that the optimal number of turns is 
aiven bv 



A' r ,„ p = \/Ri,/h T 
which yields the maximum power output 



(32) 



K 
4/1 



(33) 

Several observations can be made from the foregoing, 
viz: 

(i) The output power is proportional to of. implying 
that the specific power capability of the generator 
improves significantly as the vibration frequency 
increases. 

(ii) At the resonant frequency the plunger displacement 
Z attains its maximum value, which is proportional to 
the magnitude of the acceleration oPX of the vibration 
source. At this frequency the generator is most effective 
in converting the kinematic vibration energy into elec- 
trical energy. From this perspective, it is desirable to 
design the suspension spring such that the generator 
resonates at the frequency at which the acceleration of 
the vibration source is a maximum. 

(iii) The output power is independent of the load resist- 
ance while being dependent on K w and K, which are 
related to the electromagnetic design parameters of the 
generator. 




g- 1OO0 



1.0 o. e 

Fig. 1 0 Specific power as function of RJR, and x/R, 




0.8 



0.9 



0.6 0.7 
R m R s 

Fig. 11 Output power as function of RJR, at x f = 0.0! 1 m 

. 4.3 Design optimisation 
As pointed out earlier," the'specW'powef'capabriity of" 
the generator can be maximised by selecting a set of 
optimal design parameters such that the ratio Kj/K, 
reaches its peak for a given volumetric constraint. The 
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most influential design parameters for the machine 
topology shown in Fig. la are the ratios RJR S and 
x [ JR s and therefore an optimal design can be achieved 
by scanning the two parameters. Fig. 10a shows the 
specific output power of the generator as a function of 
the ratios RJR S and x/R„ assuming that x m h p = 0.7, 
the airgap length G = 1.0mm, Z =■ 1.4mm, R s - 
8.5mm, and / = 20 Hz, while Fig. 11 shows the genera- 
tor output power as a function of R„J£* S with x p = 
11.0mm. The specific output power reaches its maxi- 
mum at R„/R s = 0.72 and r/R, = 1.29, and is more 
sensitive to the ratio of RJR S than to rJRy 

5 Spring design 

The suspension spring of the generator fulfils two basic 
functions: to replace linear bearings that are commonly 
used in linear machines so that the plunger is free from 
mechanical friction and wear, thereby resulting in 
improved reliability and efficiency; and to facilitate 
flexible axial plunger movement with a desired spring 
constant so that the generator resonates at a specified 
frequency. These require that the spring has a very high 
radial stiffness to withstand any unbalanced magnetic 
force due to mechanical and assembly tolerances, and a 
very low axial stiffness, particularly when resonating at 
a low frequency. 

The suspension spring, shown in Fig. lb, manufac- 
tured from beryllium copper sheet using photographic 
and etching processes, can meet these requirements. 
The central hole in the spring is a precise fit onto the 
plunger shaft. The radial and axial stiffnesses of the 
spring arc affected by the thickness of the beryllium 
copper, and the number and length of the spiral 
grooves, one of whose trajectories is described by 

x = aa cos a 4- x 0 , . , OA s 
, Q-'o < « < c*i (34) 
y = aasina + y 0 

For the three-groove configuration the trajectories of 
the other two grooves are obtained by rotating eqn. 34 
by 120° and 240°, respectively. Hence, by varying the 
thickness of the sheet and the parameters a (x 0 , >*o) and 
(oq, aj) in eqn. 34, it is possible to obtain the desired 
spring characteristic, which can be predicted by nonlin- 
ear structural finite-element analysis, using the package 
Ansys for example. Fig. 12 shows the deformation of a 
spring of 0.1mm thickness, 8.5mm radius, and having 



0.32mm wide grooves, under 1.0N axial force, as pre- 
dicted by Ansys. Fig. 13 compares the predicted and 
measured spring characteristics. As seen, the predicted 
and measured results are in reasonable agreement, the 
maximum error being 27% which may be attributable 
to the fact that the material properties used in the pre- 
diction may differ from those of the actual material. It 
is also observed that the spring exhibits some nonline- 
arity, which manifests itself as an increase in the spring 
stiffness as the deflection is increased. 




0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 
axial deformation.mm 
Fig. 13 Comparison of measured and predicted spring characteristics 

predicted 

— O— measured 

Table 2: Specification and design parameters of recipro- 
cating tubular linear permanent-magnet generator 



Specification 




Design parameters 




Nominal output volt- 


5.0 


Armature radius 


8.5 


age, V 




mm 




EMF constant, V-s/m 


70.0 


Magnet radius R„, 


6.0 






mm 




Nominal deflection, mm 


*0.8 


Pole pitch t„, mm 


11.0 


Peak deflection, mm 


2.0 


Axial length of mag- 


7.5 






nets r„, mm 




Resonant frequency, Hz 


50 


Mass of plunder, g 


18.9 


Output power, mW 


20 


Airgap, mm 


1.0 



displacement, m 




displacement, m 

_ .22SE-05 

5! .S41E-03 

™ .0012" 

25 .001683 

*** .002105 



.01)2948 
.003369 
.00379 



Fig. 12 Predicted deformation of spring under 1 K axial load 
a = 0.001 ni; = .t/2: a\ = I7*/3 
a View from z-axis 
6 View from x-axis 
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Table 3: Measured and predicted winding EMF constant, resistance and inductance 



Winding parameters EMF constant V»s,/m Resistance, 20*C Inductance. H 

No. of turns Wire diam (mm) Measured Predicted Measured Predicted Measured Predicted 

5900 0.063 123 TV8 iTio 1750 0171 0.163 



6 Experimental results 

A generator, whose specification and design parameters 
arc given in Table 2, lias been prototyped, Fig. 14. Its 
winding parameters, and measured and predicted EMF 
constant, resistance and inductance arc given in 
Table 3. Measurements were carried out by mounting 
the generator to a voice-coil vibrator and applying a 
sinusoidal excitation. An accclerometer was attached to 
the vibrator to measure the input vibration. Since the 
input thrust force to the generator is proportional to 
the peak acceleration, most of the tests were conducted 
with constant peak acceleration so as to obtain com- 
parative results. 




Fig. 1 4 Prototype of reciprocating generator 



4.5 




10 20 30 40 50 60 70 

frequency f, Hz 

a 



12 i : ; — 

• 10 ; • 

cc : \ f. 

I 4 ; ;.../..;... 



10 




10 20 30 40 50 60 70 



frequency f, Hz 

Fig. 15 Open-circuit output voltage as function of vibration frequency 
uniier 0.382g peak input acceleration 

Fig. 15 shows the open-circuit output voltage as a 
function of the vibration frequency with 0.382g peak 
input acceleration. The open-circuit output voltage 
reaches its peak at a resonant frequency of 51 Hz. 
Fig. 16 shows the output voltage and power as func- 
tions of the vibration frequency with a 1.72kQ resistive 
load and l.Og peak acceleration. Comparing with 
— Fig:— HS; — afthottgh — srnrrlai — resorntnt — rx^Trriorrr'-ris" 
observed, the resonant frequency is now slightly lower 
at 46Hz. This is attributed to the influence of the non- 
linear spring characteristic shown in Fig. 13. On open- 
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Fig. 16 Output voltage and power as functions of vibration frequency 
with l.72kC2 resistive load and l.Og peak acceleration 
u Output voltage 
b Output power 

circuit, the effective damping is a minimum and a small 
input vibration may result in a larger deflection, a 
Stiffer spring behaviour and therefore a higher resonant 
frequency. Fig. 17 shows a typical output voltage 
waveform and frequency spectrum when the generator 
is vibrated at 50Hz frequency with the input accelera- 
tion waveform and frequency spectrum given- in 
Fig. 18, again with a 1.72kQ resistive load. The output 
voltage is almost sinusoidal, the harmonic distortion 
being less than 4%, which again is largely due to the 
nonlinear spring effect. While this level of distortion is 
acceptable, there arc other potential applications for 
the tubular generator topology for which the distortion 
might be one or two orders of magnitude too high. 
Velocity transducers, for example, for which the har- 
monic content could be reduced by appropriate signal 
conditioning. Fig. 19 shows the variation of the output 
voltage and power as functions of the input accelera- 
•trorrat-t h e r e so nan t fr e qu enc y of 46Hz, with ttrs -satng 
resistive load. Theoretically, the output power is pro- 
portional to the square of the input acceleration, which 
is the case when the input acceleration is lower than 
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I. Og. However, as the magnitude of the input accelera- 
tion is increased, the nonlinear stiffness of the sprin K 
restricts the attainable deformation. As a consequence, 
the subsequent increase in output power become pro- 
gressively lower than proportional to the square of 
acceleration. 




O 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 
time t,s 

b Spcctmm 



7 Conclusions 

A small reciprocating tubular linear permanent-magnet 
generator for use in applications such as on-board 
generation of electrical power on vibration transducers 
has been described, systematically analysed and 
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b Output voltage 
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experimentally characterised. A complete electrome- 
chanical model of the power generation system has 
been established, and a design methodology to achieve 
maximum output power at a specified resonant fre- 
quency has been developed. The validity and effective- 
ness of the developed analysis and design technique 
have been confirmed by finite-element analysis and 
demonstrated on a prototype generator. The techniques 
arc, by their general nature, applicable to reciprocating 
linear actuators of other topologies. 
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10 Appendix: Definition of a, n/ b| n , a„„ a„„ and B 0 

Let 

cin = B/ 0 (m-r,i?m); c 2n = BK Q (m n R m ) 
c 3n = BIi(m n R m ); c 4n = BKxim^Rra) 
c 5n = BI 0 (m n R s ); c 6 „ = BK 0 (m n R s ) 
cm = Bh(q !H R m )- c 12j = Bh{q mj R m ) 



C2u = BI 0 (qiiiR m ); c 2 2. ; = BI Q (g lnj R m ) 

C 7 n = (2/m n Tfp) [cOSTOnT// - COS ?7l n [Tj / + T m )] 

c 8n = (sin m n T U ) jm n 

Fine — — / cosm n zs\nqjn{z - T;r)dz 

7"m J 

t f + t t /2 

Fjnc ~-~r I cosm n zsmq, n j(z - r p )dz 



P^ic = — / cosqm(z - Tn)s'mm n zdz 



R nis — 



dl-n 



lip 

T. 

C27vC5n 



J cos q;iij(z - t p ) sinm n zdz 



; (kn = 



■ c 3n 



B n sin — sin — — 

mi V 2r (p J \ 2r lp 

Then a,,„ a m a Wj and B 0 are solutions of the following 
(N E + I E + J E + 1) x (N B + I E + J £ + 1) linear equa- 
tions: 

IE 

di n ai n + yXP n »C2n)qjK + c 7n B 0 

Jb 

+ ^>2(PnjsOn j )aiIIj = Bn 

N E 

53 W'ncrfan)o/n + ClUO-IIi = 0 

TV=1 

- E(^CSn)a; n + £ Sii a/ „ + ^ Bo = 0 
Nr. 

^2(Fj nc d2 n )ain +Ci2jaujj = 0 

n=l 

and 

bin = (C5n/c 67l )a;,i 
where N E) J K and J E are the numbers of the harmonic 
terms used for the calculation of the flux density in 
regions /, II and III, respectively. 
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